An autoantibody signature predictive for multiple sclerosis.

Although B cells are implicated in multiple sclerosis (MS) pathophysiology, a predictive or diagnostic autoantibody remains elusive. In this study, the Department of Defense Serum Repository (DoDSR), a cohort of over 10 million individuals, was used to generate whole-proteome autoantibody profiles of hundreds of patients with MS (PwMS) years before and subsequently after MS onset. This analysis defines a unique cluster in approximately 10% of PwMS who share an autoantibody signature against a common motif that has similarity with many human pathogens. These patients exhibit antibody reactivity years before developing MS symptoms and have higher levels of serum neurofilament light (sNfL) compared to other PwMS. Furthermore, this profile is preserved over time, providing molecular evidence for an immunologically active preclinical period years before clinical onset. This autoantibody reactivity was validated in samples from a separate incident MS cohort in both cerebrospinal fluid and serum, where it is highly specific for patients eventually diagnosed with MS. This signature is a starting point for further immunological characterization of this MS patient subset and may be clinically useful as an antigen-specific biomarker for high-risk patients with clinically or radiologically isolated neuroinflammatory syndromes.

High-throughput complement component 4 genomic sequence analysis with C4Investigator.

The complement component 4 gene loci, composed of the C4A and C4B genes and located on chromosome 6, encodes for complement component 4 (C4) proteins, a key intermediate in the classical and lectin pathways of the complement system. The complement system is an important modulator of immune system activity and is also involved in the clearance of immune complexes and cellular debris. C4A and C4B gene loci exhibit copy number variation, with each composite gene varying between 0 and 5 copies per haplotype. C4A and C4B genes also vary in size depending on the presence of the human endogenous retrovirus (HERV) in intron 9, denoted by C4(L) for long-form and C4(S) for short-form, which affects expression and is found in both C4A and C4B. Additionally, human blood group antigens Rodgers and Chido are located on the C4 protein, with the Rodger epitope generally found on C4A protein, and the Chido epitope generally found on C4B protein. C4A and C4B copy number variation has been implicated in numerous autoimmune and pathogenic diseases. Despite the central role of C4 in immune function and regulation, high-throughput genomic sequence analysis of C4A and C4B variants has been impeded by the high degree of sequence similarity and complex genetic variation exhibited by these genes. To investigate C4 variation using genomic sequencing data, we have developed a novel bioinformatic pipeline for comprehensive, high-throughput characterization of human C4A and C4B sequences from short-read sequencing data, named C4Investigator. Using paired-end targeted or whole genome sequence data as input, C4Investigator determines the overall gene copy numbers, as well as C4A, C4B, C4(Rodger), C4(Ch), C4(L), and C4(S). Additionally, C4Ivestigator reports the full overall C4A and C4B aligned sequence, enabling nucleotide level analysis. To demonstrate the utility of this workflow we have analyzed C4A and C4B variation in the 1000 Genomes Project Data set, showing that these genes are highly poly-allelic with many variants that have the potential to impact C4 protein function.

HLA-DPB1*13:01 associates with enhanced, and KIR2DS4*001 with diminished protection from developing severe COVID-19.

Extreme polymorphism of HLA and killer-cell immunoglobulin-like receptors (KIR) differentiates immune responses across individuals. Additional to T cell receptor interactions, subsets of HLA class I act as ligands for inhibitory and activating KIR, allowing natural killer (NK) cells to detect and kill infected cells. We investigated the impact of HLA and KIR polymorphism on the severity of COVID-19. High resolution HLA class I and II and KIR genotypes were determined from 403 non-hospitalized and 1575 hospitalized SARS-CoV-2 infected patients from Italy collected in 2020. We observed that possession of the activating KIR2DS4*001 allotype is associated with severe disease, requiring hospitalization (OR = 1.48, 95% CI 1.20-1.85, pc = 0.017), and this effect is greater in individuals homozygous for KIR2DS4*001 (OR = 3.74, 95% CI 1.75-9.29, pc = 0.003). We also observed the HLA class II allotype, HLA-DPB1*13:01 protects SARS-CoV-2 infected patients from severe disease (OR = 0.49, 95% CI 0.33-0.74, pc = 0.019). These association analyses were replicated using logistic regression with sex and age as covariates. Autoantibodies against IFN-α associated with COVID-19 severity were detected in 26% of 156 hospitalized patients tested. HLA-C*08:02 was more frequent in patients with IFN-α autoantibodies than those without, and KIR3DL1*01502 was only present in patients lacking IFN-α antibodies. These findings suggest that KIR and HLA polymorphism is integral in determining the clinical outcome following SARS-CoV-2 infection, by influencing the course both of innate and adaptive immunity.

Exceptional diversity of KIR and HLA class I in Egypt.

Genetically determined variation of killer cell immunoglobulin like receptors (KIR) and their HLA class I ligands affects multiple aspects of human health. Their extreme diversity is generated through complex interplay of natural selection for pathogen resistance and reproductive health, combined with demographic structure and dispersal. Despite significant importance to multiple health conditions of differential effect across populations, the nature and extent of immunogenetic diversity is under-studied for many geographic regions. Here, we describe the first high-resolution analysis of KIR and HLA class I combinatorial diversity in Northern Africa. Analysis of 125 healthy unrelated individuals from Cairo in Egypt yielded 186 KIR alleles arranged in 146 distinct centromeric and 79 distinct telomeric haplotypes. The most frequent haplotypes observed were KIR-A, encoding two inhibitory receptors specific for HLA-C, two that are specific for HLA-A and -B, and no activating receptors. Together with 141 alleles of HLA class I, 75 of which encode a KIR ligand, we identified a mean of six distinct interacting pairs of inhibitory KIR and HLA allotypes per individual. We additionally characterize 16 KIR alleles newly identified in the study population. Our findings place Egyptians as one of the most highly diverse populations worldwide, with important implications for transplant matching and studies of immune-mediated diseases.

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis.

Background: The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods. Methods: We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls. Results: We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk. Results: The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

High-throughput complement component 4 genomic sequence analysis with C4Investigator.

The complement component 4 gene locus, composed of the C4A and C4B genes and located on chromosome 6, encodes for C4 protein, a key intermediate in the classical and lectin pathways of the complement system. The complement system is an important modulator of immune system activity and is also involved in the clearance of immune complexes and cellular debris. The C4 gene locus exhibits copy number variation, with each composite gene varying between 0-5 copies per haplotype, C4 genes also vary in size depending on the presence of the HERV retrovirus in intron 9, denoted by C4(L) for long-form and C4(S) for short-form, which modulates expression and is found in both C4A and C4B. Additionally, human blood group antigens Rodgers and Chido are located on the C4 protein, with the Rodger epitope generally found on C4A protein, and the Chido epitope generally found on C4B protein. C4 copy number variation has been implicated in numerous autoimmune and pathogenic diseases. Despite the central role of C4 in immune function and regulation, high-throughput genomic sequence analysis of C4 variants has been impeded by the high degree of sequence similarity and complex genetic variation exhibited by these genes. To investigate C4 variation using genomic sequencing data, we have developed a novel bioinformatic pipeline for comprehensive, high-throughput characterization of human C4 sequence from short-read sequencing data, named C4Investigator. Using paired-end targeted or whole genome sequence data as input, C4Investigator determines gene copy number for overall C4, C4A, C4B, C4(Rodger), C4(Ch), C4(L), and C4(S), additionally, C4Ivestigator reports the full overall C4 aligned sequence, enabling nucleotide level analysis of C4. To demonstrate the utility of this workflow we have analyzed C4 variation in the 1000 Genomes Project Dataset, showing that the C4 genes are highly poly-allelic with many variants that have the potential to impact C4 protein function.

A common allele of HLA is associated with asymptomatic SARS-CoV-2 infection.

Studies have demonstrated that at least 20% of individuals infected with SARS-CoV-2 remain asymptomatic1-4. Although most global efforts have focused on severe illness in COVID-19, examining asymptomatic infection provides a unique opportunity to consider early immunological features that promote rapid viral clearance. Here, postulating that variation in the human leukocyte antigen (HLA) loci may underly processes mediating asymptomatic infection, we enrolled 29,947 individuals, for whom high-resolution HLA genotyping data were available, in a smartphone-based study designed to track COVID-19 symptoms and outcomes. Our discovery cohort (n = 1,428) comprised unvaccinated individuals who reported a positive test result for SARS-CoV-2. We tested for association of five HLA loci with disease course and identified a strong association between HLA-B*15:01 and asymptomatic infection, observed in two independent cohorts. Suggesting that this genetic association is due to pre-existing T cell immunity, we show that T cells from pre-pandemic samples from individuals carrying HLA-B*15:01 were reactive to the immunodominant SARS-CoV-2 S-derived peptide NQKLIANQF. The majority of the reactive T cells displayed a memory phenotype, were highly polyfunctional and were cross-reactive to a peptide derived from seasonal coronaviruses. The crystal structure of HLA-B*15:01-peptide complexes demonstrates that the peptides NQKLIANQF and NQKLIANAF (from OC43-CoV and HKU1-CoV) share a similar ability to be stabilized and presented by HLA-B*15:01. Finally, we show that the structural similarity of the peptides underpins T cell cross-reactivity of high-affinity public T cell receptors, providing the molecular basis for HLA-B*15:01-mediated pre-existing immunity.

Genotype List String 1.1: Extending the Genotype List String grammar for describing HLA and Killer-cell Immunoglobulin-like Receptor genotypes.

The Genotype List (GL) String grammar for reporting HLA and Killer-cell Immunoglobulin-like Receptor (KIR) genotypes in a text string was described in 2013. Since this initial description, GL Strings have been used to describe HLA and KIR genotypes for more than 40 million subjects, allowing these data to be recorded, stored and transmitted in an easily parsed, text-based format. After a decade of working with HLA and KIR data in GL String format, with advances in HLA and KIR genotyping technologies that have fostered the generation of full-gene sequence data, the need for an extension of the GL String system has become clear. Here, we introduce the new GL String delimiter "?," which addresses the need to describe ambiguity in assigning a gene sequence to gene paralogs. GL Strings that do not include a "?" delimiter continue to be interpreted as originally described. This extension represents version 1.1 of the GL String grammar.

A Predictive Autoantibody Signature in Multiple Sclerosis.

Although B cells are implicated in multiple sclerosis (MS) pathophysiology, a predictive or diagnostic autoantibody remains elusive. Here, the Department of Defense Serum Repository (DoDSR), a cohort of over 10 million individuals, was used to generate whole-proteome autoantibody profiles of hundreds of patients with MS (PwMS) years before and subsequently after MS onset. This analysis defines a unique cluster of PwMS that share an autoantibody signature against a common motif that has similarity with many human pathogens. These patients exhibit antibody reactivity years before developing MS symptoms and have higher levels of serum neurofilament light (sNfL) compared to other PwMS. Furthermore, this profile is preserved over time, providing molecular evidence for an immunologically active prodromal period years before clinical onset. This autoantibody reactivity was validated in samples from a separate incident MS cohort in both cerebrospinal fluid (CSF) and serum, where it is highly specific for patients eventually diagnosed with MS. This signature is a starting point for further immunological characterization of this MS patient subset and may be clinically useful as an antigen-specific biomarker for high-risk patients with clinically- or radiologically-isolated neuroinflammatory syndromes.

Adeno-associated virus 2 infection in children with non-A-E hepatitis.

An outbreak of acute hepatitis of unknown aetiology in children was reported in Scotland1 in April 2022 and has now been identified in 35 countries2. Several recent studies have suggested an association with human adenovirus with this outbreak, a virus not commonly associated with hepatitis. Here we report a detailed case-control investigation and find an association between adeno-associated virus 2 (AAV2) infection and host genetics in disease susceptibility. Using next-generation sequencing, PCR with reverse transcription, serology and in situ hybridization, we detected recent infection with AAV2 in plasma and liver samples in 26 out of 32 (81%) cases of hepatitis compared with 5 out of 74 (7%) of samples from unaffected individuals. Furthermore, AAV2 was detected within ballooned hepatocytes alongside a prominent T cell infiltrate in liver biopsy samples. In keeping with a CD4+ T-cell-mediated immune pathology, the human leukocyte antigen (HLA) class II HLA-DRB1*04:01 allele was identified in 25 out of 27 cases (93%) compared with a background frequency of 10 out of 64 (16%; P = 5.49 × 10-12). In summary, we report an outbreak of acute paediatric hepatitis associated with AAV2 infection (most likely acquired as a co-infection with human adenovirus that is usually required as a 'helper virus' to support AAV2 replication) and disease susceptibility related to HLA class II status.

The immunogenetics of COVID-19.

The worldwide coronavirus disease 2019 pandemic was sparked by the severe acute respiratory syndrome caused by coronavirus 2 (SARS-CoV-2) that first surfaced in December 2019 (COVID-19). The effects of COVID-19 differ substantially not just between patients individually but also between populations with different ancestries. In humans, the human leukocyte antigen (HLA) system coordinates immune regulation. Since HLA molecules are a major component of antigen-presenting pathway, they play an important role in determining susceptibility to infectious disease. It is likely that differential susceptibility to SARS-CoV-2 infection and/or disease course in COVID-19 in different individuals could be influenced by the variations in the HLA genes which are associated with various immune responses to SARS-CoV-2. A growing number of studies have identified a connection between HLA variation and diverse COVID-19 outcomes. Here, we review research investigating the impact of HLA on individual responses to SARS-CoV-2 infection and/or progression, also discussing the significance of MHC-related immunological patterns and its use in vaccine design.

Software update: Interpreting killer-cell immunoglobulin-like receptors from whole genome sequence data with PING.

Here, we demonstrate improvements to our bioinformatic pipeline, PING, which provides high-resolution genotyping of killer-cell immunoglobulin-like receptor (KIR) sequencing data, that expand the method to provide KIR interpretation from whole genome sequencing (WGS) data. We evaluated performance using synthetic sequence datasets and real-world data from the 1000 Genomes Project (1KGP). PING demonstrated high exonic genotyping performance on the synthetic sequence dataset meant to approximate real-world data at 95% accuracy (N = 1366). This result was mirrored in the analysis of 1KGP European data (N = 215) with most genes showing near or below 5% frequency of unresolved exonic genotypes, which is an important indicator for genotyping errors in real-world data. An analysis into the distributions of genotyping errors for the synthetic sequence datasets gave insights into how to further improve genotype accuracy. Similarly, an analysis into ambiguous exonic genotype frequencies for the 1KGP European data, which showed high rates of unresolved genotypes, highlighted that an effective phasing method will be an impactful future additional to the PING workflow. Together, these results demonstrate that PING can effectively provide high-resolution KIR genotyping on WGS data.

Interaction between maternal killer immunoglobulin-like receptors and offspring HLAs and susceptibility of childhood ALL.

Acute lymphoblastic leukemia (ALL) in children is associated with a distinct neonatal cytokine profile. The basis of this neonatal immune phenotype is unknown but potentially related to maternal-fetal immune receptor interactions. We conducted a case-control study of 226 case child-mother pairs and 404 control child-mother pairs to evaluate the role of interaction between HLA genotypes in the offspring and maternal killer immunoglobulin-like receptor (KIR) genotypes in the etiology of childhood ALL, while considering potential mediation by neonatal cytokines and the immune-modulating enzyme arginase-II (ARG-II). We observed different associations between offspring HLA-maternal KIR activating profiles and the risk of ALL in different predicted genetic ancestry groups. For instance, in Latino subjects who experience the highest risk of childhood leukemia, activating profiles were significantly associated with a lower risk of childhood ALL (odds ratio [OR] = 0.59; 95% confidence interval [CI], 0.49-0.71) and a higher level of ARG-II at birth (coefficient = 0.13; 95% CI, 0.04-0.22). HLA-KIR activating profiles were also associated with a lower risk of ALL in non-Latino Asians (OR = 0.63; 95% CI, 0.38-1.01), although they had a lower tumor necrosis factor-α level (coefficient = -0.27; 95% CI, -0.49 to -0.06). Among non-Latino White subjects, no significant association was observed between offspring HLA-maternal KIR interaction and ALL risk or cytokine levels. The current study reports the association between offspring HLA-maternal KIR interaction and the development of childhood ALL with variation by predicted genetic ancestry. We also observed some associations between activating profiles and immune factors related to cytokine control; however, cytokines did not demonstrate causal mediation of the activating profiles on ALL risk.

HLA variation and antigen presentation in COVID-19 and SARS-CoV-2 infection.

The extraordinary variation of the human leukocyte antigen (HLA) molecules is critical for diversifying antigen presentation to T cells. Coupled with the rise of novel strains and rapidly evolving immune evasion by SARS-CoV-2 proteins, HLA-mediated immunity in COVID-19 is critically important but far from being fully understood. A growing number of studies have found the association of HLA variants with different COVID-19 outcomes and that HLA genotypes associate with differential immune responses against SARS-CoV-2. Prediction studies have shown that mutations in multiple viral strains, most concentrated in the Spike protein, affect the affinity between these mutant peptides and HLA molecules. Understanding the impact of this variation on T-cell responses is critical for comprehending the immunogenic mechanisms in both natural immunity and vaccine development.

Minimal observed impact of HLA genotype on hospitalization and severity of SARS-CoV-2 infection.

HLA is a critical component of the viral antigen presentation pathway. We investigated the relationship between the severity of SARS-CoV-2 disease and HLA type in 3235 individuals with confirmed SARS-CoV-2 infection. We found only the DPB1 locus to be associated with the binary outcome of whether an individual developed any COVID-19 symptoms. The number of peptides predicted to bind to an HLA allele had no significant relationship with disease severity both when stratifying individuals by ancestry or age and in a pooled analysis. Overall, at the population level, we found HLA type is significantly less predictive of COVID-19 disease severity than certain demographic factors and clinical comorbidities.

A short HLA-DRA isoform binds the HLA-DR2 heterodimer on the outer domain of the peptide-binding site.

The human leukocyte antigen (HLA) locus encodes a large group of proteins governing adaptive and innate immune responses. Among them, HLA class II proteins form α/ß heterodimers on the membrane of professional antigen-presenting cells (APCs), where they display both, self and pathogen-derived exogenous antigens to CD4+ T lymphocytes. We have previously shown that a shorter HLA-DRA isoform (sHLA-DRA) lacking 25 amino acids can be presented onto the cell membrane via binding to canonical HLA-DR2 heterodimers. Here, we employed atomistic molecular dynamics simulations to decipher the binding position of sHLA-DRA and its structural impact on functional regions of the HLA-DR2 molecule. We show that a loop region exposed only in the short isoform (residues R69 to G83) is responsible for binding to the outer domain of the HLA-DR2 peptide-binding site, and experimentally validated the critical role of F76 in mediating such interaction. Additionally, sHLA-DRA allosterically modifies the peptide-binding pocket conformation. In summary, this study unravels key molecular mechanisms underlying sHLA-DRA function, providing important insights into the role of full-length proteins in structural modulation of HLA class II receptors.

Multiple sclerosis therapies differentially affect SARS-CoV-2 vaccine-induced antibody and T cell immunity and function.

BACKGROUNDVaccine-elicited adaptive immunity is a prerequisite for control of SARS-CoV-2 infection. Multiple sclerosis (MS) disease-modifying therapies (DMTs) differentially target humoral and cellular immunity. A comprehensive comparison of the effects of MS DMTs on SARS-CoV-2 vaccine-specific immunity is needed, including quantitative and functional B and T cell responses.METHODSSpike-specific Ab and T cell responses were measured before and following SARS-CoV-2 vaccination in a cohort of 80 study participants, including healthy controls and patients with MS in 6 DMT groups: untreated and treated with glatiramer acetate (GA), dimethyl fumarate (DMF), natalizumab (NTZ), sphingosine-1-phosphate (S1P) receptor modulators, and anti-CD20 mAbs. Anti-spike-Ab responses were assessed by Luminex assay, VirScan, and pseudovirus neutralization. Spike-specific CD4+ and CD8+ T cell responses were characterized by activation-induced marker and cytokine expression and tetramer.RESULTSAnti-spike IgG levels were similar between healthy control participants and patients with untreated MS and those receiving GA, DMF, or NTZ but were reduced in anti-CD20 mAb- and S1P-treated patients. Anti-spike seropositivity in anti-CD20 mAb-treated patients was correlated with CD19+ B cell levels and inversely correlated with cumulative treatment duration. Spike epitope reactivity and pseudovirus neutralization were reduced in anti-CD20 mAb- and S1P-treated patients. Spike-specific CD4+ and CD8+ T cell reactivity remained robust across all groups, except in S1P-treated patients, in whom postvaccine CD4+ T cell responses were attenuated.CONCLUSIONThese findings from a large cohort of patients with MS exposed to a wide spectrum of MS immunotherapies have important implications for treatment-specific COVID-19 clinical guidelines.FUNDINGNIH grants 1K08NS107619, K08NS096117, R01AI159260, R01NS092835, R01AI131624, and R21NS108159; NMSS grants TA-1903-33713 and RG1701-26628; Westridge Foundation; Chan Zuckerberg Biohub; Maisin Foundation.

Severe delayed hypersensitivity reactions to IL-1 and IL-6 inhibitors link to common HLA-DRB1*15 alleles.

OBJECTIVES: Drug reaction with eosinophilia and systemic symptoms (DRESS) is a severe, delayed hypersensitivity reaction (DHR). We observed DRESS to inhibitors of interleukin 1 (IL-1) or IL-6 in a small group of patients with Still's disease with atypical lung disease. We sought to characterise features of patients with Still's disease with DRESS compared with drug-tolerant Still's controls. We analysed human leucocyte antigen (HLA) alleles for association to inhibitor-related DHR, including in a small Kawasaki disease (KD) cohort. METHODS: In a case/control study, we collected a multicentre series of patients with Still's disease with features of inhibitor-related DRESS (n=66) and drug-tolerant Still's controls (n=65). We retrospectively analysed clinical data from all Still's subjects and typed 94/131 for HLA. European Still's-DRESS cases were ancestry matched to International Childhood Arthritis Genetics Consortium paediatric Still's cases (n=550) and compared for HLA allele frequencies. HLA association also was analysed using Still's-DRESS cases (n=64) compared with drug-tolerant Still's controls (n=30). KD subjects (n=19) were similarly studied. RESULTS: Still's-DRESS features included eosinophilia (89%), AST-ALT elevation (75%) and non-evanescent rash (95%; 88% involving face). Macrophage activation syndrome during treatment was frequent in Still's-DRESS (64%) versus drug-tolerant Still's (3%; p=1.2×10-14). We found striking enrichment for HLA-DRB1*15 haplotypes in Still's-DRESS cases versus INCHARGE Still's controls (p=7.5×10-13) and versus self-identified, ancestry-matched Still's controls (p=6.3×10-10). In the KD cohort, DRB1*15:01 was present only in those with suspected anakinra reactions. CONCLUSIONS: DRESS-type reactions occur among patients treated with IL-1/IL-6 inhibitors and strongly associate with common HLA-DRB1*15 haplotypes. Consideration of preprescription HLA typing and vigilance for serious reactions to these drugs are warranted.

A common allele of HLA mediates asymptomatic SARS-CoV-2 infection.

Despite some inconsistent reporting of symptoms, studies have demonstrated that at least 20% of individuals infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) will remain asymptomatic. Although most global efforts have focused on understanding factors underlying severe illness in COVID-19 (coronavirus disease of 2019), the examination of asymptomatic infection provides a unique opportunity to consider early disease and immunologic features promoting rapid viral clearance. Owing to its critical role in the immune response, we postulated that variation in the human leukocyte antigen (HLA) loci may underly processes mediating asymptomatic infection. We enrolled 29,947 individuals registered in the National Marrow Donor Program for whom high-resolution HLA genotyping data were available in the UCSF Citizen Science smartphone-based study designed to track COVID-19 symptoms and outcomes. Our discovery cohort (n=1428) was comprised of unvaccinated, self-identified subjects who reported a positive test result for SARS-CoV-2. We tested for association of five HLA loci (HLA-A, -B, -C, -DRB1, -DQB1) with disease course and identified a strong association of HLA-B*15:01 with asymptomatic infection, and reproduced this association in two independent cohorts. Suggesting that this genetic association is due to pre-existing T-cell immunity, we show that T cells from pre-pandemic individuals carrying HLA-B*15:01 were reactive to the immunodominant SARS-CoV-2 S-derived peptide NQKLIANQF, and 100% of the reactive cells displayed memory phenotype. Finally, we characterize the protein structure of HLA-B*15:01-peptide complexes, demonstrating that the NQKLIANQF peptide from SARS-CoV-2, and the highly homologous NQKLIANAF from seasonal coronaviruses OC43-CoV and HKU1-CoV, share similar ability to be stabilized and presented by HLA-B*15:01, providing the molecular basis for T-cell cross-reactivity and HLA-B*15:01-mediated pre-existing immunity.